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Abstract 

Using a rational i?-matrix associated with the 4x4 defining matrix represen- 
tation of C2 = sp(4), the Lie algebra of 5^(4), a one-site operator solution of the 
associated Yang-Baxter algebra acting in the Fock space of two harmonic oscillators 
is derived. This is used to define iV-site integrable systems, which are soluble by a 
version of the algebraic Bethe ansatz method without nesting. All essential aspects 
JC I of the work generalise directly from C2 to Cn ■ 

o' 
o. 

5 1 Introduction 

G^ , One can see from references such as [I|] that the algebraic Bethe ansatz continues to be 

a subject of interesting development and improvement. This paper is devoted to applying 
the method to integrable models with underlying Lie algebras other than oi = su{2) in a 
way that produces unnested Bethe equations. 
O ■ Let V and Ti denote auxiliary and quantum spaces respectively. Define R{u) with 

spectral parameter u m.V ®V, and T{u) in V C?) 7i. Then the Yang-Baxter equation for 
R{u) on V ® V V, 



i?12(n - t;)i?i3(n)i?i3(^) = i?23Mi?13 (n)i?12(n " v) , (1) 

and the Yang-Baxter operator algebra relation on V ® V ® ?^, 

Ri2{u - v)T^{u)T2iv) = T2iv)T^{u)Ruiu - v) , (2) 

are objects of central importance to the study of integrable systems. In contexts in which 
there is an underlying algebra of type a„, especially Oi, things are quite well understood. 
See e.g. [^, |^, Q. However, even for a„ for n > 1, solution of related models by the 
Bethe ansatz method is iterative, giving eigenstates specified by means of nested Bethe 
equations 0. The same applies also to models with so{n) or c„ = sp{2n) as underlying 
Lie algebra [0, |^. A paper that describes the nested Bethe ansatz method clearly is 
0. Papers which state the case for seeking, and which implement, alternatives to that 
method with its evident complications are |]I0| (on Gaudin models) and ||11|| . 



Here we deal with theories governed by c„ = sp{2n). We first obtain a solution of 
(|l|), and then a new 'one-site' solution L{u) = T{u) of (Q) in V ® ?i when 7i is the 
representation space of the metaplectic representation M.^ of c„ defined in terms of 2n 
pairs of harmonic oscillator variables. This is the main result of the paper. Given this 



L{u) however, we can proceed directly down the path laid down by Sklyanin ^ for the 
construction of families of iV-site integrable models. 

The integrable models just indicated are of interest because they are soluble by means 
of a one stage application of the Bethe ansatz method, giving rise to Bethe equations 
without any nesting, equations that are no more complicated than those familiar within 
oi studies p. 

In this paper, we present some essential results for C2 and indicate their extension to 
c„ in section 2, describe A^„ in section 3, and derive our formulas for R{u) and L{u) in 
section 4. Section five contains a sketch of the solution of our C2 integrable models by our 
Bethe ansatz method. Fuller discussion of this last topic and various related matters will 



be provided in papers in preparation |12 



2 The Lie algebra Cn 

We present the basis of C2 in Cartan-Weyl form with simple roots ri = (1, — 1), r2 = (0, 2), 
so that the other positive roots are ra = ri + r2 = (1, 1) and r4 = ri + rs = (2,0). Let 
H = {Hi, H2) denote the generators of the Cartan subalgebra of C2. Let -E±a, a G {1, .., 4}, 
denote the raising and lowering generators. Then the Lie algebra C2 is specified by 

[Ei,E2] = V2Es , [EuEs] = V2E^, etc. (3) 

We also use a Cartesian basis Xi,i G {1, . . . , 10}, for C2 

Hi = Xi , H2 = X2 \'2E±a = X2a+1 + iX2a+2 ■ (4) 

Then we may define the structure constants of C2 via 

[Xi,Xj\ =icijkXk , (5) 

and the quadratic Casimir is 

C'^'^=X,Xi = H^i+Hi + Y,{Ea,E^^} . (6) 

a 

We also require the matrices of the four- by-four defining representation T> of C2, defined 
according to 

Xi^ Xi ; Hi ^ hi , H2 ^ h2, E±a ^ e±„ . (7) 

We can display these all at once by use of a 4 x 4 matrix C = Xi® Xi = XiXi that is very 
useful for work on the C2 Yang-Baxter algebra. We have 

4 
C = Hihi + H2h2 + Y,iEae-a + E_^e^) 

/ Hi E_i E_3 72^-4 \ 

El H2 v2-E'-2 E^3 

E3 V2E2 —H2 —E^i 

V V2E, E; -El -Hi J 



One may read the explicit forms for hi,h2, e±a off (@); verify that they do represent (^) 
correctly, and note the hermiticity properties h\ = hi, hi = hi, e^„ = e^a appropriate to 
generators of the compact representation X> of C2. It can further be seen that the matrices 
Xi,i G {1,2. ..10}, possess the properties 

xl = Xi, Tr (xj) = 0, Tr (xiXj) = 26ij, xj = -JxiJ^ . (9) 

Here J, the standard C2 symplectic form, in our representation, has, as its only non-zero 
elements: Ju = J23 = —J32 = —J^i = i- In view of (j^), it follows that the matrices Xi 
satisfy the completeness relation 

XiacXibd = SadSbc " JabJcd Or Xi®Xi = P -K , (10) 

where Kab,cd = JabJcd, is proportional to the projector onto the symplectic trace, and P 
is the permutation map on C^ ^ C^. 

In the defining representation Xi \-^ Xi the Casimir C*-^-* has eigenvalue 5. In the 
adjoint representation (Xj)_,fc = —icijk, it has eigenvalue 12, so that also 

CijkCiji = I25ki . (11) 

Our work on the Yang-Baxter algebra requires us to extend the span of the Xi to the space 
of all traceless hermitian 4x4 matrices. Thus we define a specific set of five matrices ya, 
a G {1, 2, ..5}, with the properties 

yi = ya, Tr?/a = 0, TTyayb = 25ab, Va = JVaJ'^ , (12) 

so that the Jxi define a set of 10 symmetric matrices, whereas Jya and J itself define a 
set of 6 antisymmetric ones. Further we can write 

XiXj ~r XjXi Oij ~r ^(^ijaVa , v J 

which defines d^a = djia such that dua = 0. Use of the completeness relation (|^) and 
trace properties, allows proof of the results 

(^ija^^ijb "J^ab ; Va "^(^ijaXiXj . \^^) 

Our definition ( p^ of the i/a refiects the fact that the set A^ = {xi, ...,Xio,yi, ...,1/5} can 
serve as a set of 15 a^ Gell-Mann A-matrices which satisfy 

X\ = \a, TrAA = 0, Ti Xa\b = 26ab, A^ ® A^ = 2P - I/4 , (15) 

where I4 is the 4x4 unit matrix. If we write X = Xi®Xi,Y = ya^ya, then the last part 
of (|T3p and ( PUJ ) read as 

X + Y = 2P-lh, X = P-K , (16) 

which can be solved for P and K in terms of 14 , X , Y. 

3 The Metaplectic Representation, A4n of Cn. 

Just as the metaplectic representation of ci = ai is discussed using the Fock space of one 
harmonic oscillator ||13|, n so also is that of c„ discussed using n independent oscillators. 



We present results explicitly for C2 in a form that generalises completely and naturally for 
c„ for all n. 

Since we define the row vector f ^ in terms of two sets of standard harmonic oscillator 
variables by 

v^ = {a[ al 02 ai) , (17) 

then, using also the definition of J given above, it can be shown explicitly that the 
operators 

X, = \v'^x,Jv (18) 

obey the same commutation relations (|^) as the matrices Xi of (^. Then (^) leads us to 
the explicit results 

Hi = -{ai, a\}, H2 = |{a2, 4} , E+i = a\a2, E^i = ai4 , 



E+2 = -al^/V2, E+3 = -a\al E+^ = -af/V2, (19) 

E_2 = o^ll^i -^-3 = ai«2, E^4^ = a\/\/2 . 



One can also check directly that ( [l9l) correctly provides a representation of (^. From 
(p!9|), we find the hermiticity conditions H{ = Hi, Hi = H2, E\_-^ = -E-i, E'l = 
—E^a, a G {2,3,4}. These reflect the fact that ([T9|) generates the infinite dimensional 
unitary metaplectic representation A^2 of the real non-compact groups Sp{2,M.). Also 
the maximal compact subalgebra of (|19D is that of an SU{2) x f/(l) group, with angular 
momentum type generators given by 

J^ = \{Hi-H2), J+ = a{a2 = E+i = A, J^=a\ai = E^i (20) 

and f/(l) generator A = \{Hi + H2), [A, J] = 0. In addition, we note the two commut- 
ing ai subalgebras generated by Hi,Ej-4 and H2,E±2- These correspond to the two ci 
metaplectic subrepresentations of our C2 metaplectic representation ([T9|). 

We also need to verify that the 'non-compact' raising and lowering generators of C2 
can be arranged to define vector operators ^ and x^ using the standard Racah definition 
of tensor operators. Explicitly ^ and x have spherical components 

(e-i,eo,e+i) = (^2,^3,^4) , (x-i,xo,x+i) = (^-4,-^-3,^-2) . (21) 

All the required properties can be checked explicitly. It is significant for the Bethe ansatz 
solution of the integrable models we construct, that the components of the vectors ^ 
commute with each other, as do those of x- 

We are now in a position to derive the key algebraic result for 7^2, a result which 
allows us to obtain an infinite dimensional operator solution of (|l|) when A^2 is taken as 
the quantum space 7i, (and similiarly for c„). To this end we use (^B|) and ( pUD to write 
the matrix C = x-iXi of (^) in the form 

C = -\h + {Jv)v^ , (22) 

whence it follows that C obeys the quadratic relation 

C^ = -3C - \h . (23) 



We stress that this non-trivial result holds for the particular (metaplectic) representation 
(P^) of C2- It does not hold in an arbitrary representation of C2. 

Since TrC = 0, (|2|) implies that for M2 we have C^^^ = ^TrC^ = -|, a result that 
can be checked directly to be correct by inserting (|TU|) into (^), and expressing everything 
in terms of Ni = a|ai, N2 = a2a2- 

The entire discussion generalises readily from C2 to c„. The generalisation of (|1^), and 



of the completeness relations ([T^) are obvious, and imply that the only change needed in 
(p^) replaces the unit matrix 1 4 of C^ by l2n- Then we have 

C = -{n + 1)C - {'-^) hn , (24) 

so that 

C(2) = lTrC2 = -idim(c„). (25) 

4 The Yang-Baxter Equation and Algebra of Cn 

Let V = C^. Then we seek a solution of the Yang-Baxter equation (|1|) of C2 of the form 
R{u) =uI + r]P + C{u)K, or 

RabM"^) ='^^aJbd + V^adScb + C{u)JabJcd ■ (26) 

This is a natural generalisation of the work for a2n where only I and P are invariant under 
the action of sl{2n, C), to the present case where K also is invariant under the action of 
sp{2n, C). 

Equation (^) yields a solution of (|1]) provided that C{u) is given by 

C{u) = -^^, X = -3v , (27) 

A — u 

and the same holds for c„ when A = —{l + n)ri. These results are the symplectic analogues 



of results found for the orthogonal groups, e.g. in |T5[, and agree with results presented 
in different form, e.g. in 0. We wish to use the i?-matrix (|26|) , ( ^Tf ) to reach solutions 
of the Yang-Baxter algebra on V ® V ® 7Y when 7i is the module corresponding to the 



operator representation given by (19 



First as a guide we review the route followed successfully and in generality for the 
same purpose in the case of the a„ algebras. There one has 



R{u) = uI + r]P = uI + 7]Uj2^A®XA + ^l] , 



(28) 



using the completeness relation of the Gell-Mann A-matrices of a„. Changing the repre- 
sentation ^\a to Xa leads to 

L{u) = {uI + ^7])+vJ2^^^^^ ■ (29) 

A 

which can be seen to give a 'one-site' solution T{u) = L{u) of (^. 

The same approach does not work straightforwardly for C2 (or c„) because (]26|) involves 
not only X = Xi® Xi but also Y = ya® Va, where (c/. section 2) the i/a do not belong to 



the Lie algebra C2- However, as we have shown above, they can be expressed in terms of 
the Xi and so we enforce the change of representation from Xi to Xi onto Y 

ya®ya=ya® \dijaXiXj ^ Va® \dijaXiXj (30) 

= i(-5,, + {x,,x,})®X,X, = i(C2 + 3C-iC(2)) . 

where, as before, C = Xi ® Xi = XiXi , C^^^ = XiXi. The resuhs (0), (0) and (|) 
have all been used here. We do not know how to handle the complication that follows the 
general use of (pOD, but we do know how to proceed in special cases, such as that of section 
3, in which C*^^^ = — | and (p3D holds, for then the RHS of (0) vanishes. Accordingly we 
seek a solution T{u) or first L{u) of (0) in the form 



This ansatz succeeds for 



L{u) = a{u)I + l3{u)xiXi , (31) 

a{u) u — 5 



/3{u) T] 



6 eR , (32) 



and we may put (3{u) = rj and a{u) = u — S. We emphasise the non-trivial nature of 
the result obtained. Since R{u) given by ( ^61) necessarily involves Y = ya®ya the ansatz 
(pl|) is not obviously valid a priori in our work. It seems unlikely to be valid for general 
representations in the quantum space. But it works for metaplectic representation of C2, 
and likewise for other c„. 

5 Bethe Ansatz for C2 Integrable Models 

Sklyanin's general procedure |Q| for constructing integrable spin models, e.g. if operators 
Xi = Si acting in ?-^ = V^ where V^ is the vector space in which the spin-s representation 
of SU{2) acts, is available here also once we know a basic ('one-site') matrix L{u) G V^Ti 
which satisfies (^. We write 

N 

T{u)=ICY[L{u~6r) , (33) 

where /C is a constant matrix such that [/C® /C, R{u)] = 0. Since T{u) is an n-site solution 
of (^, it follows in a well-known way that [t{u),t{v)] = where t{u) = Tr T{u), and 
expansion of t{u) in powers of u yields constants of the motion. 

In this section we sketch the solution for the eigenstates of t{u) = '^i^i Tu of the 
system described by (^) when L{u) is the one-site solution of (|^) given in the last section 
by means of the algebraic Bethe ansatz method 0, ^. Although detailed attention to 
/C and to the Sr is important for the discussion of the completeness of the set of these 
eigenstates, we defer this to a later publication, here setting /C = 1 and 6r = 0. We 
wish rather to explain how the Bethe ansatz method works in our models, emphasising 
the simplifications that stem from special features of the representation ( [T9| ) of Xi, and 
exploiting the transformation properties of the Tij under the compact su{2) subgroup of 
C2. The latter follow for the Tij when J = Ylr=i "^r is used, J^ for r = 1, . . . A^ being 
defined at the r-th site in terms of the oscillator variables of that site as in (|19D . Perhaps 
it should be pointed out that J as just defined is not equal to the (A^ — l)-fold coproduct 
of the 'one-site' operator (7+ = T21 , Jz = KT^ii ~ ^22 ) ' ^- ~ ^12 ); i^ the sense of 



the Yang-Baxter algebra @. Nevertheless its commutation relations with the A'^^^T^- , 
involving the standard comultiplication, can be shown explicitly to give these operators 
the correct su{2) tensor operator properties. 

Our sketch deals only with the tower of Bethe states based in the Fock vacuum of 
both the oscillators at each of the A^ sites. Similar work for other towers whose ground 
states involve oscillator states of occupation number one requires only a modest extension 
of the analysis described here. See [0. We note first, from (p!9D, that all the Lij with 



i < j, as well as L21 and L43 annihilate the Fock vacuum at each site, and it is easy to 
promote the same result to the same set of A^-site Tij of (pS]). Also, at any site we have 
(L42 — L3i)|0) = 0, and this result to can be similarly promoted to 

(T42-T3i)(m)|0) = . (34) 

We aim first to construct the Bethe states 

Tiiiv)\0) , T42(t;i)T3i(t;2)|0) , ... , (35) 

which have j = m = 1,2, ... quantum numbers w.r.t J. Then the other states with the 
same j and lower m follow by application of the lowering operator J_. Since [J, t{u)] = 
all the states of each multiplet of Bethe states have the same eigenvalue of t{u), and 
allowed f-values given by the same set of Bethe equations. One naturally expects (and 
finds) that there are, alongside the j = 2 multiplet, j = 1 and j = multiplets bilinear in 
the creation operators T4i,T3i,T42 and T32. There is not scope in this paper to describe 



in full the subtleties of the analysis (in which (|3^) allows significant simplifications), so 
next we present results. 

Let w = u — v,\et t+{u) denote the eigenvalue of Tn and of T22 for the vacuum state 
|0), and let t_(m) do the same for T33 and T44. It is easy to calculate t±{u) using 
and (|l9|) . Then for j = 1, we find 



t{u)T,iiv)\0) = r{u)T,i{v)\0) , (36) 

where the eigenvalue t{u) is given by 

r(«) = (l + ^K(«) + (l-^K(«) , (37) 

and V is determined by the Bethe equations t+{v) = r_(f) 

For j = m = 2, the eigenvalue t{u) of t{u) for r4i(t>i)T4i(f2)|0), and the Bethe 
equations for vi and f 2 are given by 



r{u) = 2(l + ^ + ^ + 2^^)r+(n) 

+ 2(1 - JL ^ JL + 2^^)t4u) , (38) 



where we write Wk = u — Vk^ k = 1^2 and v = vi — f2, and 

T+{vi)/r_{vi) = T_{v2)/r+{v2) = {v + 2r])/{v - 2r]) . (39) 

The generalisation to higher integral j- values is evident. Turning to m = 1 states bilinear 
in the creation Tij, we find two, one that follows by application of J_ to the state just 
discussed. Writing 



T42(t;i)T4i(t;2)|0) , 

T3i(t;i)T4i(t;2)|0) , 

T,iiv,)T,2iv2)\0) , 

niivi)T3iiv2)\0) , (40) 



we find the |21) state is given by {(f>i + </>2 + 03 + (p^), and the unique |11) Bethe state 
in the context is given by {(0i + 02 — 03 — 04) + ^(0i ~ 02)} with t{u) eigenvalue given 
by 2(1 + ^ + :^)t+{u) + 2(1 - ^ - ;^)r_(M), and Vi and V2 determined by the Bethe 
equations t+{vi) = T_(fi) ,T+(t>2) = r_(t>2). The pattern of the Bethe states should by 
now be clear. 

The states constructed have various nice properties. Their Bethe equations ensure 
that the poles of their t{u) eigenvalues have residues zero. They are also invariant under 
the exchange of vi and V2, although proof of this is not obvious, requiring detailed use of 
the commutation relations of the Tij that we are using. 

A brief sketch of how these results are derived is called for, because there is a level 
of complication not seen in Bethe ansatz studies in models with oi = 5^(2) invariance. 
This is no doubt the price to be paid for obtaining results without recourse to a nesting 
process. Complications not present in ai studies come into play here whenever the third 
term of ( p6D enters the commutation relations amongst the Tij non-trivially. A sufficient 
illustration of what has to be done whenever this happens is provided by showing how 
to compute the effect of the term Tii{u) of t{u) on T4i(t>)|0). We here shall abbreviate 
any Tij{u) or Tij{v) by T^j or T[-. Our strategy is push to the right end of any term only 
factors that annihilate |0), or else are of the type Tkk-, with no sum on fc, which multiply 
|0) by T±. Direct calculation of wTuT'^]^ from (|^) fails to achieve this, leaving a term 
in T2iT^i that has to be eliminated by calculation from (|^) of wT2iT^i. This two stage 
procedure yields finally the result 

TnT^, = (l + ^)T;iTn-— T4ir{i + ^(T^iT2i-r3ir2'i) . (41) 

W WW 

Another example yields 

TuTi, = (1 + ^)ri,Tn + (^ - ^)riiri2 - ^T4it;, 

w + r] w w + T] w 

- ^T,,T^, + {^-^)T',,T22 + ^n,T2, . (42) 

w w w + r] w + 7] 

One must push each of the four pieces of t{u) past each of the creator Tij, needing a total 
of about 32 results in all, to deduce results such as those just quoted. This entails eight 
results like each of (^) and (|42|); the remainder do not involve third term of ( pG] ) and are 



written down directly. 

6 Discussion 

We have described, for C2, formalism and the construction of a family of integrable models 
that generalise in every respect to all c„. This apphes to the construction (pTS]) and the 
proof that the Xi obey the same commutation relations as the Xi, to the deduction of 
relations ( |3TD and (p^). A similar construction is available for the dn series, where we 
use fermionic instead of bosonic oscillator variables. Note that in the c„ case, there is 
a maximal compact u{n) subalgebra, and sets of creation operators Tij which transform 
according to its dimension n defining representation (1, ... ,0), leading to a u{n) multiplet 
structure of Bethe eigenstates |T2[ , with a pattern similar to that indicated in section five. 



We note also that our method of getting Bethe eigenstates without recourse to nesting 
applies equally well to the treatment ||12|] of c„ Gaudin models |l^ . 
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